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ABSTRACT 

Radio recombination lines (RRLs) are powerful, extinction-free diagnostics 
of the ionized gas in young, star-forming regions. Unfortunately, these lines 
are difficult to detect in external galaxies. We present the results of EVLA 
observations of the RRL and radio continuum emission at 33 GHz from NGC 
253, a nearby nuclear starburst galaxy. We detect the previously unobserved 
H58« and H59a RRLs and make simultaneous sensitive measurements of the 
continuum. We measure integrated line fluxes of 44.3 ± 0.7 W mr 2 and 39.9 ± 
0.8 W mr 2 for the H58a and H59a lines, respectively. The thermal gas in NGC 
253 is kinematically complex with multiple velocity components. We constrain 
the density of the thermal gas to 1.4 — 4 x 10 4 cm -3 and estimate an ionizing 
photon flux of 1 x 10 53 s -1 . We use the RRL kinematics and the derived ionizing 
photon flux to show that the nuclear region of NGC 253 is not gravitationally 
bound, which is consistent with the outflow of gas inferred from the X-ray and Ha 
measurements. The line profiles, fluxes, and kinematics of the H58a and H59a 
lines agree with those of RRLs at different frequencies confirming the accuracy of 
the previous, more difficult, high frequency observations. We find that the EVLA 
is an order of magnitude more efficient for extragalactic RRL observations than 
the VLA. These observations demonstrate both the power of the EVLA and the 
future potential of extragalactic RRL studies with the EVLA. 

Subject headings: Galaxies: ISM — Galaxies: star formation — radio continuum: 
galaxies — radio lines: galaxies — Galaxies: individual (NGC 253) 
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Introduction 



How stars form is one of the most important physical questions in astrophysics today. 
While we have viable theories that explain how s tar formation occurs in the Milky Way 
( IZinnecker fc Yorkd 120071 ; iMcKee fc Ostrikerl 120071 ). we do not yet understand how stars 
form in galaxies with different conditions than those in the Milky Way. A key piece of the 
puzzle is quantitative information about the physical conditions in the interstellar medium 
(ISM) in extragalactic star-forming regions with a wide range of metallicities, total masses, 
and star formation rates. Studies of young star-forming regions still embedded in their natal 
clouds of dust and dense gas are particularly important; these embedded regions give us 
crucial information about the ISM before the newly-formed stars have destroyed their birth 
environments and thus tell us about the conditions in which stars form. However, these 
regions are poorly understood because they are heavily obscured in the optical and near- 
infrared (Ay > 10), which limits the utility of "standard" optical/IR nebular diagnostic 
lines. 

Radio recombination lines (RRLs) are powerful tracers of the physical properties of 
obscured thermal gas. These li nes occur when electrons make transitions between high levels 
(n > 40) of an excited atom (IRoelfsema fc Gosd Il992l ; iGordon fc Sorochenkd 120091 ) . For 
typical extragalactic spatial res olutions, RRLs a ct like density filters with higher frequency 
lines tracing higher density gas ( iZhao et al.lll997l ). Observations of RRLs at multiple, widely 
spaced frequencies constrain the density and filling factor of the thermal gas in a star- 
forming region without being affected by the significant extinction from the dust and gas 
within the star-forming region. RRLs also provide kinematic information about the thermal 
gas and estimates of the ionizing photon flux, the star formation rate, and the star formation 
efficiency Unfortunately, the use of RRLs outside the Milky Way system has been limited 
to a dozen bright, nearby galaxies because these lines are faint (a few percent above the 
continuum) and have large velocity widths (> 100 kms -1 ). Moreover, the higher frequency 
RRLs, which trace the higher density thermal gas and thus younger star forming regions, 
have been particularly challenging to observe because the bandwidths needed to observed 
these lines were larger than the maximum available bandwidth of the previous generation of 
instrument at ion . 

Work by Anantharamaiah and collaborators established the V ery Large Array (VLA) 



as a powerful instrument for extragalactic RRL observations fe.g.. lAnantharamaiah et al. 



1993 



Zhao et al 



1996 



Anantharamaiah fc GosJ 1996; Zhao et al. Ii99t1: 



Phookim et al 



1998; 



Anantharamaiah et al. 



2000 



iDas et al. 



2001 



Mohan et al.l2001 



2002 



et all2005: Mohan et al 
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20051 : IRov etaD \200M Most of these observations were at frequencies < 8.3 GHz; this was 
the highest frequency for which the VLA correlator had a large enough velocity coverage 
to observe the entire line using a single correlator tuning. Observations of RRLs at higher 
frequencies (43 GHz) by Rodriguez- Rico and collaborators required combining four corre- 
lator tunings in post-pr ocessing for the necessary velocity coverage ([Rodriguez- Rico et al. 
20041 . bood . I2006l . l2007al lbh . While those observations provided vital information on the high 
density gas in several systems, they may introduce systematic uncertainties in the line profile 
and in the continuum subtraction, particularly for very broad lines. 

The bandwidth o f the National Rad io Astronomy Observatory (NRAO) Expanded Very 
Large Array (EVLA; Perley et al. 2011 1^1 WIDAR correlator and the improved sensitivity 
and additional frequency coverage of the EVLA receivers are opening new frontiers in RRL 
studies. These improvements will make observations of these lines significantly more efficient 
than with the VLA enabling us to broaden the sample of galaxies with RRL measurements. 
We present observations taken as part of the EVLA Resident Shared Risk Observing (RSRO) 
program AK726 to test the RRL capabilities of the EVLA. The electronics of the EVLA are 
completely new, so the performance of the EVLA for detecting these lines cannot be directly 
extrapolated from the very successful attempts using the VLA. 

For our pilot study, we chose NGC 253, a nearby galaxy with a nuclear st arburst 
(jUlvestad fc Antonuccil 119971 ; I Johnson et al.l l200ll ; iFernandez-Ontiveros et al.l 120091 ) . This 
galaxy is an ideal source to test the RRL capabilities of the EVLA because of its strong, 



well-studied KR 



(|Seaauist fc Bell 



1997 



emission; it has measurements of 11 different RRLs from HI 66a to H40a 
Mebold et alil98dlAnantharamaiah fc Gos^ll99oiE996llPuxlev etal. 



1977 



Mohan et al. 



2002 



2005 



Rodriguez- Rico et al.ll2006l ). We present observations of two, 



previously undetected RRLs (H58a and H59a). We demonstrate that the EVLA produces 
results comparable to previous VLA observations in an eighth of the time and allows simulta- 
neous observations of multiple RRLs. With the bandwidth of the EVLA WIDAR correlator, 
we are able to more accurately derive the line profile and make excellent, simultaneous mea- 
surements of the continuum. We use the H58a and H59ck lines to constrain the density of 
the thermal gas in NGC 253, estimate the ionizing photon flux, and explore the kinematics 
of the thermal gas in the nuclear region. 

We assume a distance to NGC 253 of 3.4 Mpc (1" is 16 pc), which is the mean of the 
distances determined from three different fields using the tip of the red giant branch method 



1 Roy et al. ( 20051 2008 ) also include data from the Australia Telescope Compact Array. 

2 The National Radio Astronomy Observatory is a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. 
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flDalcanton et alJl200a l 



2. Data Reduction 

The data were taken with the EVLA in DnC configuration on 2010 September 22, 25-30. 
The angular resolution was approximately 1.8" with a field of view of 1.4'. The resolution 
corresponds to a physical resolution of 30 pc, which is five to ten times the size of a typical 
massive star cluster. While these observations cannot resolve individual clusters, they are 
also unlikely to resolve out diffuse emission associated with the clusters. 

We observed H58a (32.85220 GHz) and H59a (31.22332 GHz), which were previously 
undetected in NGC 253. We configured the correlator to provide two basebands, each con- 
taining eight adjacent 128 MHz wide sub-bands with 128 channels and dual polarization. 
We tuned the basebands to the sky frequencies corresponding to H58a and H59a for NGC 
253. The sensitivity at the edges of the sub-bands decreases, which increases the noise in 
this region of the spectrum by a factor of ~ 1.4. The placement of the line peak at the edge 
of a sub-band results in the noise being higher in the region with the line peak. However, 
this placement allows very accurate measurements of the flux in the line wings. 

We used J0137+3309 (3C48) and J0542+4951 (3C147) as absolute flux density calibra- 
tors, J0319+4130 as the bandpass calibrator, and J0137-2430 as the complex gain calibrator. 
The total on-source time was 57 minutes. 

We used CASA§| to calibrate the data using the standard EVLA high frequency calibra- 
tion methocEl with three modifications. First, since at the time these data were calibrated 
there were no Ka-band flux density calibrator models, we used the K-band flux density 
calibrator models. The next two modifications were used to ensure that each individual 
sub-band was on the same flux density scale. The average opacity value for the central 
two spectral windows was used as the opacity value for all spectral windows in a baseband. 
Since the opacity values typically vary by only 2-3% across the entire baseband, this does 
not introduce any significant systematic effects. Finally, since the flux density derived for the 
complex gain calibrator was the same (within the errors) for each sub-band in one baseband, 
we set the flux density of the complex gain calibrator for each sub-band to be the average 
of the flux densities determined for the complex gain calibrator. 



http : / / casa . nrao . edu 
4 http : / / casaguides . nrao . edu 
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We self-calibrated using the continuum emission deriving one solution for each baseband. 
We continuum subtracted the data in u — v space and combined the line data (shifting to 
the same reference frequency) and the continuum data. Table [T] lists the image properties. 



Results 



Accurate measurements of both the RRL and continuum emission are necessary to model 
the RRL emission. Previous RRL observations sacrificed continuum sensitivity to detect and 
spectrally resolve the RRL. The bandwidth of the EVLA allows us, for the first time, to make 
sensitive, simultaneous measurements of both the RRL and continuum emission. Figure [T] 
shows the 33 GHz continuum and integrated H58a line flux density for NGC 253 overlaid on 
an HST image of the optical and IR emission. We do not show the 32 GHz or the H59a data 
since they are very similar. The strongest RRL emission is coincident with the peak of the 
radio continuum emission. There is anothe r extended region of RRL em ission toward the 
southwestern portion of the disk. Following IRodrfguez-Rico et al.l ( 120061 ). we refer to these 
regions as the northeastern (NE) and southwestern (SW) regions. We measured the flux in 
the NE region using a rectangle with a bottom left corner of (00 h 47 m 33^3, — 25°17'19) and a 
top right corner of (00 h 47 m 33^0, — 25°17'15) and the SW region using a rectangle with a bot- 
tom left corner of (00 h 47 m 33!l, -25°17'23) and a top right corner of (00 h 47 m 32!8, -25°17'19). 

The extinction in this region is extremely high. In the HST image, blue (B-band) 
represents unobscured young star formation, while the red (Paschen a) traces obscured 
recent star formation. The RRL emission is coincident with the highly extincted recent 
star formation south of the dust lane. The higher resolution of the HST data allows us to 
distinguish many of the individual star formation regions that are not resolved in the RRL 
observations. Even at the larger spatial resolution of the EVLA data, we can distinguish 
two separate, obscured star-forming regions associated with primary and secondary radio 
continuum and RRL peaks. 

The distribution of the H58a line and 33 GHz continuum emission are similar to the 



distrib ution of the H53a line and 43 GHz continuum determined by IRodrfguez-Rico et al. 



(120061 ) . Both data sets have similar noise levels. The observing times required to reach these 
sensitivities, however, are drastica l ly dif ferent. The on-source observing time for the data 
presented in iRodriguez-Rico et al.l (120061 ) was 8 hours, while the H58a d ata was obtained in 
57 minutes. The observations presented in IRodrfguez-Rico et al.l (120061 ) only measured one 
line and a small amount of continuum, while we measured two lines as well as obtaining ex- 
cellent continuum data at both line frequencies. This comparison shows that the EVLA is at 
least an order of magnitude more efficient than the VLA for extragalactic RRL observations. 
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The 8 GHz bandwidth of the full correlator will deliver even higher observing efficiencies by 
providing simultaneous observations of at least four RRLs at these frequencies. 



Figure [2] shows channel maps for the H58a emission. The NE and SW regions have 
different kinematics. The NE region covers a wide range of velocities from ~ kms -1 to 
360 kms^ 1 to, while the SW region covers a much narrower range in velocities (190 kms 
to 300 kms " 1 ). The central velocities of the NE and SW region are offset by approximately 



60 kms . lAnantharamaiah fc Gossl ( 19961 ) attribute the complex H92a velocity field to 



either a tumbling bar potenti al, the merger or accr e tion o f another galaxy, or a warp in 
the disk in the nuclear region. iRodriguez-Rico et al.l (120061 ) suggest that the RRL velocity 
structure could also be the res ult of a outflow driven by the nuclear starburst. T here is addi- 
tional X-ray and Ha evidence ([Strickland et al.ll2000l ; IWestmoquette et al.ll201ll ) supporting 
a nuclear outflow. 

Figure [3] shows the spectrum for the entire line emitting region and for the NE and 
SW regions. As expected from the models of these lines, the higher frequency (H58a) 
line has a higher peak flux density than the lower frequency line (H59a). The H58a and 
H59c>! line profiles for the e ntire RRL region are similar to the line profiles reported by 
Rodriguez-Rico et al.l (120061) an d show the presence of several kinematic components dis- 
cussed in detail by lMohan et al.l (120051 ). We averaged every three channels in the H58a cube 
together to compare the peak of the H58a line profile to the peak of the H53a line profile from 
Rodriguez-Rico et al.l (120061 ). The line profiles agree well. The H53a peak line flux density 
is larger than the H58ct peak line flux density, again in agreement with phenomenological 
expectations. 

We fit Gaussians to the H58a and H59a line profiles. While Gaussians are not the best 
fit for the data because of the multiple velocity components, these fits allow us to directly 
compare our results to the fits for other RRLs in th e literature. Table [2] gives the parameters 
for the fit to the lines as well as the fits presented in lRodriguez-Rico et al.l (120061 ). The H58a; 
and H59a measurements in this paper agree with the measurements of the line flux from 
the literature. Differences between the flux measurements for each region in the current and 
archival data can be attributed to differences in the regions selected. The differing velocity 
widths of the H58a and H59a are due to systematic uncertainties in the Gaussian fits because 
of multiple velocity components and the different sensitivities of the two images. The signal 
to noise of the H58a data is higher than that for the H59a data because the H59a line is 
fainter than the H58a line and the noise in the H59a cube is slightly higher. 

Extragalactic RRL emission is modeled as a collection of HII regions, each with an identi- 
cal electron density ( n P ), size (7), and temperature (T), immersed in a region of diffuse radio 



continuum emission (jPuxley et al.l 11991 



Anantharamaiah et al.l Il993l . 120001 ) . We use the 
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method of IZhao et al.l (119971 ) to estimate the range of densities consistent with the H58a line 
data to compare with those derived in the literature using other RRLs. Assuming a typical 
non-thermal spectral index (-0.7) and a temperature of 10 4 K, H58a line flux and continuum 
flux density are consistent with electron densities of 1.4 — 4 x 10 4 c m~ 3 . This corresponds 



to an ionizing photon flux of 10 x 10 52 s l . Using Equation (1) in lAnantharamaiah et al. 



(120001 ). we obtain star formation rates averaged over the lifetime of OB stars of 2 M© yr 1 in 
the RRL emission region. The electron density, ionizing photon fluxes, and s tar formation 



rate e stimates derived from the H 58a line agree with those in the literature (IMohan et al. 



20051 ; iRodriguez-Rico et al.ll2006l ). Our ionizing photon flux is slightly 



than the 3 — 6 x 10 52 s 1 ionizing photon flux derived by IMohan et al.l (120051 ) 



ligher 10 x 10 5 s 



We use the derived ionizing photon flux and the measure d line widths to es timate 
whether the nucleus of NGC 253 is gravitationally bound (e.g. Turner et al. l2003h . We 



estima te the mass of the cluster from the number of ionizing photons using iLeitherer et al. 
(199 9J). We assume that the clu ster was formed in an instantaneous burst, is 3 Myr old 



(IFernandez-Ontiveros et al.l |2009| ). and has a Salpeter IMF with an upper mass limit of 



100 Mfn, which gives a clust er mass of 4.0 x 10 6 M . The mass of the ionized gas derived by 
Rodriguez-Rico et al.l (120061 ) is negligible compared to the cluster mass. The line emitting 
region is 53 pc by 135 pc. The escape velocity for the nucleus is then between 17 kms -1 and 
26 kms -1 . This width is much smaller than the velocity width of the RRL emission. We 
assume that the RRL line widths are dominated by the motion of the gas since they are much 
larger than the intrinsic width of a RRL from a single HII region (~ 20 kms" 1 ). The escape 
velocity is still much smaller the RRL velocity width if we consider only the NE region. 
There we derive an ionizing photon flux of 9.3 x 10 52 s _1 , a cluster mass of 3.97 x 10 6 M and 
an escape velocity of 25 kms -1 . These calculations show that the central region of NGC 253 
is most likely not gravitationally bound, supporting the hypothesis that there is an outflow 
of gas from the nucleus of this galaxy. 

Although RRLs provide an extinction-free tracer of the gas ionized by the cluster, the 
use of RRL data may underestimate the number of ionizing photons if the photons are 
absorbed by dust rather than ionize hydrogen. Therefore, the mass of the cluster is a lower 
limit. To obtain an escape speed comparable to the RRL velocity width, the number of 
ionizing photons would have to be underestimated by a factor of 45. 



4. Conclusions 



We have presented observations of the H58a and H59a RRLs and the continuum emis- 
sion at 32 and 33 GHz from the nucleus of NGC 253. We measure integrated line fluxes of 
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44.3 ±0.7 and 39.9 ±0.8 W m -2 for the H58a and H59a lines respectively. We constrain the 
thermal gas density to 1.4 — 4 x 10 4 cm -3 . We derive an ionizing photon flux of 10 x 10 52 s" 1 . 
These estimat es of the density and ionizing photon flux are consist ent with estimates from 
the literature jMohan et alil2002[ I2OO5I : iRodriguez-Rico et al-lhood ). 



The kinematics of the thermal gas in NGC 253 are complex. The NE region has a 
broader velocity width (170 - 200 kms" 1 ) than the SW region (145 - 159 kms" 1 ). From 
our estimate of the ionizing photon flux, we estimate a cluster mass for the nuclear region of 
NGC 253 of 4 x 10 6 M . The escape velocity for this mass is only 17 kms -1 to 26 kms -1 : 
much lower than the H58a and H59a velocity widths. Therefore, the nuclear region of NGC 
253 is most likely not gravitationally bound. We suggest that an outflow of gas from the 
nuclear region is the most likely explanation for the therm al gas kinematics of N GC 253. 
The presence of this outflo w has been inferred from X-ray ( {Strickland et al.ll2000f ) and Ha 
(IWestmoquette et al.ll201ll ) observations. 



In contrast to previous VLA high frequency RRL observations, we were able to detect 
not one, but two RRLs, and make much more sensitive measurements of the continuum in 
an eighth of the time on source. Accurate measurements of both the RRL and continuum 
emission are necessary to model the RRL emission from these regions in detail. The EVLA 
is currently an order of magnitude more efficient for high frequency RRL observations than 
the VLA. The EVLA observing efficiency will increase even further once the full 8 GHz 
bandwidth is available. 

These observations demonstrate the potential of the EVLA to revolutionize extragalactic 
RRL observations. Future EVLA observations will be able to use RRLs to measure the 
properties of the thermal gas in a wider variety of galaxies. Ultimately, these observations 
will lead to a quantitative picture of how the properties of the ISM affect the resulting stellar 
populations. 



The authors acknowledge the hard work and dedication of the EVLA commissioning 
team. K.E.J, acknowledges support provided by NSF through CAREER award 0548103 and 
the David and Lucile Packard Foundation through a Packard Fellowship. 
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Table 1. Image Properties 
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Table 2. RRL Parameters 







Wrest 


Sc 


S L 


AVfwhm 


Vhel 


Sl AVfwhm 




Region 


Line 


GHz 


mjy 


mjy 


km s _1 


km s — 1 


10" 22 Wm- 2 
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Fig. 1.— HST ACS (F475W and F814W) and NICMOS (F187N) images of NGC 253. The 
F187N filter (Paschen a) is red, the F814W filter (I band) is green, and the F475W filter 
(B band) is blue. Left: 33 GHz continuum contours from 5 to 500 a on a logarithmic scale 
(la = 0.09 mJy beam -1 ). Right: Integrated H58a emission contours on a logarithmic scale 
between 5% and 90% of the peak flux of 1.15 Jy beam" 1 kms" 1 . The distribution of the 
32 GHz continuum and H59a line emission is identical. The beam (~ 1.8") is shown in the 
lower left hand corner of each panel. 
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Fig. 2. — Channel images for H58a emission. For this figure, we have averaged every three 
channels together. The contours are -3a, 3a, 5a, 7a, 9a, etc (la = 0.43mJy beam -1 ). The 
upper right hand corner of each panel gives the heliocentric velocity (optical definition). The 
lower left hand corner of the first panel shows the beam (~ 1.8"). 
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Fig. 3. — The H58a (left panels) and H59a (right panels) emission for the entire RRL region 
(top row), the NE region (middle row), and the SW region (bottom row). The positions 
of each region are indicated in Figure [TJ The horizontal line indicates a flux density of 
mJy. The dotted vertical line indicates the boundary between two sub-bands. The line in 
the upper right hand corner of each panel is the 3a per channel error for each region . The 
dashed line in the top row is the H53a line profile from iRodrfguez-Rico et al.l ( 120061 ). The 
velocities are heliocentric (optical definition). 



